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Abstract 

Double-exchange model interacting with adiabatic phonons is studied in the presence of randomness by using 
the Monte-Carlo method and the systematic size-scaling analysis. A bicritical behavior is found between the 
ferromagnetic metal and the charge-ordered insulator. We find the distinct response to the randomness between 
these two states, which agrees well with the experimental results in the colossal magnetoresistance manganites. 
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Phase competition and randomness play a key 
role in the colossal magnetoresistance (CMR) 
phenomena in manganese oxides [1,2]. The for- 
mer phase competition between the ferromagnetic 
metal (FM) and the charge-ordered insulator 
(CO) leads to large fluctuations near the multicrit- 
ical point [3]. The latter randomness suppresses 
the transition temperatures above which CMR 
is much enhanced [4]. Theoretically, except for a 
phenomenological argument [5], this multicritical 
phenomenon has not been fully investigated thus 
far. Therefore, it is strongly desired to study a 
realistic model for this phenomenon including the 
metal-insulator transition. 

In this work, we study the multicritical phenom- 
ena in the double-exchange (DE) model interact- 
ing with phonons whose Hamiltonian reads 



(1) 



Here, the first line is the DE part which is derived 
in the limit of the large Hund coupling with the 
Ising symmetry (Si = ±1) [6], and the second line 
represents the electron-phonon coupling and the 
elastic energy of phonons. We consider the classi- 
cal phonon in the breathing mode. The A term is 
for the cooperative distortion of the MnOg tetra- 
hedra. The first term in the last line denotes the 
random on-site potential. We consider here the 
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Fig. 1. The Binder parameter and the numerical derivatives 
of the temperature dependence of the order parameter, (a) 
and (b) are for FM at g = 0.8, and (c) and (d) are for CO 
at g = 1.2, respectively. Critical temperatures are indicated 
by the crossing of the Binder parameter in (a) and (c), and 
the peak of the numerical derivatives in (b) and (d). 

binary-type distribution Si ~ ±A. The last term 
is a Gintzburg-Landau-type term which describes 
the competition between FM of the DE origin and 
CO of the checker-board type, whose order param- 
eters are denoted as M and X, respectively. We 
consider that this term mimics the elements be- 
yond this simple model such as the orbital degen- 
eracy and the antiferromagnetic coupling between 
localized spins. We consider model (1) at the elec- 
tron density n = 0.5 on the square lattice. We set 
the halfband width 14^ = 4t = 1 as an energy unit, 
and take A = 0.1 and 7 = 1 in the following. 

We study the thermodynamics in model (1) by 
MC method [7] . We have typically 10000 samplings 
for measurement after 1000 steps for thermaliza- 
tion in the absence of the randomness. In the pres- 
ence of the randomness, the measurement is typ- 
ically performed 1000 times for a configuration 
of the random potential {si}, and the results are 
random-averaged for typically 16 different config- 
urations. The transition temperatures for FM and 
CO are determined by the Binder parameter as 
well as the inflection point of the temperature de- 
pendence of each order parameter. Figure 1 shows 
a part of these analyses (A = 0). Both results con- 
sistently indicate the emergence of the long-range 
ordering. 

Figure 2 displays the phase diagram of model 
(1). In the pure case A = 0, the phase diagram 
shows the bicritical behavior where FM and CO 
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Fig. 2. Phase diagram of model (1). 

phases touch with each other at g ~ 0.9. When we 
introduce the randomness, these two phases show 
a distinct response; FM is rather robust whereas 
CO is surprisingly fragile to the randomness. Con- 
sequently, the CO phase is destroyed rapidly, and 
instead the FM state extends to larger-5 regions. 
This is a consequence of the phase competition 
near the bicritical point; if one phase diminishes, 
the other takes over. These aspects well agrees 
with the experimental results in the CMR man- 
ganites [4]. Therefore, we beheve that model (1) 
captures the essential physics of this multicritical 
phenomenon in manganites. In experiments, the 
enhanced CMR is observed near the region where 
FM replaces CO [4]. Transport property and elec- 
tronic structure in the present model will provide 
a key to understand the mechanism of CMR. 
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